Photon echo quantum memories in inhomogeneously broadened two level atoms 
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Here we propose a solid-state quantum memory that does not require spectral holeburning, instead 
using strong rephasing pulses like traditional photon echo techniques. The memory uses external 
broadening fields to reduce the optical depth and so switch ofi' the collective atom-light interaction 
when desired. The proposed memory should allow operation with reasonable efficiency in a much 
broader range of material systems, for instance Er^+ doped crystals which have a transition at 
1.5 ^m. We present analytic theory supported by numerical calculations and initial experiments. 
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I. INTRODUCTION 

Photon echo based techniques, particularly in rare 
earth ion dopants, have been long investigated for classi- 
cal signal processing More recently they have found 
use in quantum memory applications and are now leading 
the field with the longest storage times Jl, Q , the high- 
est efficiencies 0| , broadest bandwidths |5|| and the high- 
est time-bandwidth products 6']. The photon echo tech- 
niques used for quantum information applications differ 
from those used in the classical domain. Classical infor- 
mation processing using rare earth ion dopants relies on 
strong optical pulses to cause rephasing of the atomic co- 
herences. Applying a similar approach in the quantum 
realm would be highly desirable, but analysis has shown 
that the basic techniques, namely the two-pulse and 
three-pulse |9i] photon echo sequences, will not work as 
quantum memories. 

This work proposes low noise photon echo techniques 
based on optical rephasing pulses that could be useful 
as quantum memories. We present analytical theory and 
initial experiments. 

The two approaches currently used for photon echo 
quantum memories techniques, controlled reversible in- 
homogeneous broadening (CRIB) 0, [TMl and atomic 
frequency combs (AFC)[il 13 

both require spectral 
holeburning. This is the process where the inhomoge- 
neously broadened absorption profile is modified by fre- 
quency selective optical pumping to shelving states. 

For CRIB, spectral holeburning is used to create a 
spectrally narrow absorption feature. This is then broad- 
ened with a field gradient to accept the input light. The 
inverse spectral-width of the initial feature determines 
the storage time, and how far it is broadened determines 
the bandwidth. The light is then recalled by reversing the 
field gradient. To do this with high efficiency, significant 
optical thickness is required for the broadened feature. 
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which means that extremely large optical thicknesses are 
required for the unbroadened feature. Hedges et al.[l] 
started with an absorption of -130 dB in the initial fea- 
ture to achieve 69% efficiency with a time-bandwidth 
product only large enough to faithfully store one pulse. 

In AFC protocols the material is prepared in a similar 
way to CRIB except that a number of narrow absorption 
features are created. Using this technique large band- 
width delays m an d reasonable efficiencies have also been 
demonstrated 17]. 

With both CRIB and AFC methods, high efficiencies 
require efficient optical pumping and this requires long- 
lived shelving states. To date this has meant using shelv- 
ing states in the dopant ions' hyperfine structure. Our 
proposal does not rely on spectral holeburning, and the 
removal of this restriction is very appealing. The hy- 
perfine splittings are generally less than 100 MHz, which 
makes large bandwidth operation for the AFC and CRIB 
protocols problematic. Sagiamyurek et al. [s] used a 
trick to overcome the problem of hyperfine structure and 
bandwidth with AFC memories. By choosing the cor- 
rect comb spacing, they had the anti-holes that result 
from spectral holeburning appearing on top of the comb 
teeth. The downside of this approach is that high finesse 
frequency combs are not possible, severely limiting the 
efficiency. 

Another reason to avoid spectral holeburning is the 
requirement of efficient holeburning limits the choice 
of material systems. For example, while single photon 
memories have been demonstrated at telecommunication 
wavelengths using erbium [l^, the lack of an efficient 
holeburning mechanism severely limited the memory ef- 
ficiency. 

Ensemble quantum memories rely on the fact that the 
collective interaction with an optical field is greatly en- 
hanced for systems with reasonable optical depths. The 
problem with optically rephased memories is that the 
strong optical rephasing pulses necessarily invert the pop- 
ulation, this causes gain and the associated noise. 

Here we propose avoiding this problem by using a per- 
turbing field to spectrally broaden the ensemble while in 
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the troublesome collective interaction between the light 
and the atoms is effectively turned off while the atoms 
are in the excited state. 

While the dephasing caused by the applied broadening 
stops the two pulse echo forming, so long as the broaden- 
ing applied before and after the first 7r-pulse are balanced 
an echo of the 2PE will form after a second vr-pulse. This 
is the output from the memory. We call this type of echo, 
which involves both static and controlled inhomogeneous 
broadening a hybrid photon echo (HPE). 



III. ANALYTICAL THEORY 



FIG. 1: (Color online) Sequence of pulses used to create a 
hybrid echo. The dashed line shows where the 2PE would 
form if there was no electric field applied to the sample. Also 
shown are the three time regions used when we explain the 
echo sequence. 



the excited state to the point where it is optically thin. 
This removes the collective enhancement and stops the 
unwanted noise processes. We call this method where 
both reversible inhomogeneous broadening and optical 
rephasing methods are used hybrid photon echo rephas- 
ing (HYPER). 

Another approach that avoids the use of spectral hole- 
burning by using Raman echoes has also been proposed 



II. DESCRIPTION 

The hybrid photon echo rephasing technique uses a 
combination of broadening field pulses and optical pulses 
to rephase a small optical input pulse. For what follows 
we shall assume that the broadening field is an electric 
field gradient which causes inhomogeneous broadening 
due to the atom's linear Stark effect. Using a linear elec- 
tric field gradient is not a fundamental requirement of 
the HYPER protocol. As long as the broadening is com- 
pletely reversed any applied field will be acceptable. The 
case of a linear field gradient is considered here as this is 
what was used in the experimental work. 

Using two TT-pulses allows the ions to rephase near the 
ground state, forming an echo of the normal two pulse 
echo (2PE). This is beneficial as an echo forming in the 
excited state is inherently noisy. The electric field gradi- 
ents are applied at appropriate times to eliminate excited 
state noise on the output. 

The pulse sequence, shown in Fig. [1] has three time 
regions. In the first region (0 < i < ^2) one or more 
input pulses are applied followed by the application of 
a broadening field. A tt pulse is applied at ^2- The tt- 
pulse creates gain which a 2PE would ordinarily propa- 
gate through and thus undergo amplification. The effect 
of broadening the ensemble is to reduce the gain, and the 
two pulse echo, to an arbitrarily low level. In doing so 



Here we use the same theoretical description that was 
developed for rephased amplified spontaneous emission 
[7:]. It is assumed that the input light to all regions is 
weak compared to a 7r-pulse and that the rephasing pulses 
are ideal 7r-pulses. 

Figure [T] shows the pulse sequence with three different 
time regions. The equations of motion for time region 1 
before the Stark shift is applied are as follows, 
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Here, a{D) is the quantum optical (atomic) field op- 
erator, subscript 1 indicates region 1, a is the optical 
depth parameter and A represents the detuning of the 
atom. The solution for the atomic field (Eq. Hal) is given 
by 
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To obtain the solution for the optical field we Fourier 
transform Eq. [2]with respect to time and substitute into 
the transformed version of Eq. Ilbl The solution in the 
time domain is found to be 
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where 'Di{z,t) is defined as ^ Di{z,t = 0, A = u) 

We now have the complete solution in region 1 for before 
the Stark shifting fields are applied. 

We now assume that the Stark shifting field is turned 
on just before the first 7r-pulse and is strong and tem- 
porally short, allowing the dynamics of the optical field 
to be ignored. The Stark shifting field alters an atom's 
detuning dependent on its position by ijz such that 
A — A + i]z, whereas the 7r-pulse inverts the atoms lead- 
ing to Z) — >■ D'l'. The Stark shift and 7r-pulse result in 
the following transformation on the atomic field at the 
region 2 {t ~t2) boundary, 
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where rji = ri{tf—ti). Here, rj is the intensity of the Stark 
shifting field and tf — ti is the duration of that field. 

Under the small input pulse approximation the atoms 
in region 2 are all near the excited state. 

A temporally long Stark shifting pulse is applied dur- 
ing all of region 2. The equations of motion are given by 
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where again the subscript denotes region. It is noted 
that these equations are similar in form to the region 1 
equations, with the atomic fields daggered due to the tt- 
pulse and the detuning is shifted with a field intensity 
•q' such that A A + 77'z. Also, a sign change in front 
of the optical field in the atomic equation of motion is 
present. 

The atomic field solution for region 2 is given by 
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Using the region 2 boundary condition stated earlier, the 
optical field solution for region 2 is found in a similar 
fashion as for region 1: 
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This completes the solution for region 2. It is noted 
that when the Stark fields are set to zero, the two 
pulse photon echo solutions are retained with efficiency 
??2P_E = 4sinh^(^) as expected. 

We can now form the region 3 boundary condition, 
in a similar fashion to the region 2 boundary condition, 

namely 1)3(2, i = t4, A) = D\{z^ti^A) . 

The equations of motion that describe the dynamics 
of region 3 are exactly those stated for region 1, Eqs. [Tal 
and llbl Hence the optical solution for region 3 is identical 
in form to Eq. |3l with subscript 1^3 and T)3{z,t) = 

J-^ [1)3(2, ^4, a;)e~''^**j . For balanced Stark fields, the 



output optical solution in region 3 is 
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where we have defined the operator 

63(0, i) = 03(0, *) + ia / dz'e'^f>i{z\t-2t4 + 2t2). 
Jo 

In the limit of large Stark shift, Eq. [S] becomes 
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The solution has three terms. The first term is the in- 
evitable optical vacuum input at the region 3 time bound- 
ary which decays as a function of the optical depth. The 
second term is the HYPER echo which forms at a time 

= ti -I- 2t4 - 2*2 with ^HPE = (az)2e-"^ This echo 
has a maximum efficiency of 54% at a/ = 2. The third 
term contains the atomic degrees of freedom. If counter- 
propagating TT-pulses are used, causing the echo to be 
emitted in the backward direction, 100% efficiency is pos- 
sible (see appendix). 

It can be seen that the Stark shifting field over region 
2 eliminates contributions from region 2 on the region 3 
output. Taking the limit of infinite Stark shift is phys- 
ically equivalent to decoupling the optical and atomic 
fields in the equations of motion in region 2, thus 'switch- 
ing off' the collective atom-light interaction in this region 
and reducing the noise on the output field. In the decou- 
pled regime, the atomic boundary condition at region 3 
(t = ^4) is found by propagating the region 2 boundary 
condition {t = 12) forward in time by — t2. The result 
obtained for the output field in region 3 is identical to 
Eq.H 



IV. ECHO DEMONSTRATION 

The experimental setup used is shown in Fig. [21 a mod- 
ified Coherent 699 dye laser with sub 5 kHz linewidth 
drives the ■^H4 — ^D2 transition, at 605.977 nm, in a 
Pj.3+.Y2Si05 crystal. The sample is oriented such that 
the laser propagates along the crystal b-axis, and is lin- 
early polarized along the D2-axis. 

The crystal is surrounded by four electrodes in a 
quadrupole arrangement that enabled an electric field 
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FIG. 2: (Color online) Experimental setup. PBS - Polariz- 
ing beamsplitter, NPBS - non-polarizing beamsplitter, AOM - 
acousto optic modulator. A third AOM was used (not shown) 
before the interferometer to provide the required frequency 
shifts for repumping the praseodymium dopants to the re- 
quired hyperfine states. The 0.02% Pr^+:Y2Si05 crystal is 
mounted between four electrodes in a quadrupolar arrange- 
ment that are used to apply the electric field gradients to the 
samples. 



gradient to be produced along the direction of propa- 
gation. The electric fields used in region 1 varied from 
±15 V/cm over the optical path length, and were di- 
rected parallel to the b-axis. Using the tabled value of 
the Stark shift [20], this results in detunings of approx- 
imately ±2 MHz. The direction of the applied elec- 
tric field was parallel to the crystals b-axis. The entire 
sample/electrode setup is mounted inside a cryostat and 
cooled to 3. IK. Two acousto-optic modulators (AOMs) 
were used to gate the laser beam to produce the input 
pulses. The AOMs also introduce a net 10 MHz fre- 
quency shift, enabling the light exiting the crystal to be 
detected with phase sensitive heterodyne detection. 

Figure [3] shows the intensity of the light exiting the 
sample during a hybrid photon echo sequence. Three op- 
tical pulses are applied, the input pulse and two 7r-pulses. 
The input pulse is a Gaussian with a full width at half 
maximum (FWHM) power of 1.8/is, and the transmit- 
ted portion appears at (i). In the absence of the Stark 
field, the first 7r-pulse (ii) causes a 2PE (in). The 2PE 
is rephased to form the hybrid photon echo at (vi). The 
two TT-pulses were not placed symmetrically about the 
2PE so that the three pulse photon echo (3PE) at (v), 
due to the input and two 7r-pulses, is separated in time 
from the HPE. 

The amplitude of the pulses were calibrated by apply- 
ing a 2PE sequence with the first pulse half the length of 
the second. The lengths that maximised the size of the 
echo were taken as being a 7r/2 pulse and tt pulse. 

Applying the electric field to the sample causes the in- 
tensity of the 2PE to be reduced by > 99% through the 
mechanism explained earlier. The 3PE also experiences a 
large reduction in intensity, which can be explained using 
similar reasoning. The HPE is still formed albeit with re- 
duced efficiency. We believe this reduction in efficiency is 
most likely due to imperfect balancing of the broadening 
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FIG. 3: (Color online) An experimental hybrid photon echo 
trace (a) without and (b) with the electric field gradient ap- 
plied. The two TT-pulses saturated the detector. 
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FIG. 4: (Color online) Disrupting the two pulse echo with the 
broadening field. The energy in the two pulse echo is plotted 
versus the magnitude of the electric field gradient applied to 
the sample. 



applied before and after the first 7r-pulse. 

Figure 2] shows the intensity of the 2PE as the size 
of the broadening field is increased. The electric field 
varies linearly along the propagation axis, resulting in a 
top-hat distribution of frequency shifts. Thus the phase 
shifts of the ions will have a top-hat distribution as the 
2PE forms. This leads to a sinc-squared behaviour in the 
energy of the 2PE, analogous to a single slit diffraction 
pattern. 

Due to optical pumping distributing population 
amongst the hyperfine levels Pr^+:Y2Si05 exhibits spec- 
tral holeburning [21,] . Optical repumping was used so 
that the number of ions and thus the optical depth re- 
mained consistent between shots. The laser was first 
scanned over the spectral region of interest to burn a wide 
spectral hole. This hole was then filled using laser light 
pulses detuned by a specific combination of the hyper- 
fine splittings. Varying the number of these repumping 
pulses allowed for variation of the optical depth. 

The efficiency of the 2PE and HPE was measured over 
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FIG. 5: (Color online) Efficiency of the 2PE and HPE as the 
optical depth is varied, (a) Zero applied electric field, and (b) 
9 V applied to electrodes. 
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FIG. 6: (Color online) Area of the hybrid photon echo as 
the area of the input pulse is varied. 9V was applied to the 
electrodes. 



a range of optical depths (see Fig. ^ . The maximum ef- 
ficiency for the 2PE occurs at an optical depth of around 
1.5, while the maximum efficiency of the HPE is at an op- 
tical depth of around 1.2. This is independent of whether 
the electric field was on or off. The maximum 2PE effi- 
ciency reduces from 40% to less than 0.3% when 9 V was 
applied to the electrodes. 

In order to verify the linearity and determine the ef- 
ficiency of the HPE, the area of the input pulse was 
compared to the area of the output pulse (see Fig. 
The length of the input pulse was set to 1.8/is (FWHM), 
whilst the amplitude of the pulse was altered by adjust- 
ing the radio frequency power applied to A0M2. The 
optical depth for this experiment was set to 1. As shown 
in Fig. ini the relationship between input pulse and out- 
put pulse is linear for input pulses small compared to a 
TT-pulse and then flattens off for more intense pulses. 



NOISE MEASUREMENTS 



To characterize the noise properties, the heterodyne 
signal of the HYPER pulse sequence was recorded with 
no input pulse, both with and without the electric field 
applied. The Pr'^+:Y2Si05 sample used for the noise 
measurements was 4x4x20 mm (20 mm along the beam 
propagation direction), doped at 0.005%. Before the ini- 
tial TT-pulse, 100 V is applied to the two electrodes. In 
this alternate sample geometry this results in an elec- 
tric field varying from -30 V/cm to -1-30 V/cm along the 
20 mm axis, corresponding to a frequency shift of ap- 
proximately ±4 MHz. This voltage is decreased to 37.5 V 
for the longer period between the 7r-pulses to maintain a 
zero net Stark shift such that an input pulse would be 
rephased as a HPE. 

The noise was calculated as follows. For each of 
~30 000 shots, the time-domain signal was multiplied by 
the temporal envelope of the would-be 2PE (FWHM of 
1.8 ^s, centered 15 /xs after the first 7r-pulse), numerically 
heterodyned (optical beat frequency here was 6 MHz) 
down to DC and summed. This quantifies the field ampli- 
tude in the 2PE mode for each shot. Taking the variance 
across all shots yields a value related to the average field 
intensity in this specific mode. For normalization pur- 
poses, the shot-noise (noise level with local oscillator on, 
but signal beam blocked) was measured as the variance 
in an equivalent mode, but with the temporal envelope 
shifted to before the first 7r-pulse. 

At an optical depth of 0.15, the quadrature variances 
in the echo mode without the electric field applied were 
measured as 8.64±8% normalized to the shot-noise level. 
With the electric field applied, the noise reduces approx- 
imately two orders of magnitude to 1.045±4%, a level in- 
distinguishable from the vacuum. Figure [7] shows the av- 
eraged time dependence of the noise in the spectral mode 
of the echo. The experiment was repeated for larger opti- 
cal depths of 0.36 and 0.98, where the noise was similarly 
reduced by factors of 150 and 160 respectively. 

There are three mechanisms through which the applied 
field gradient decreases the noise. Firstly, the gain fea- 
ture created by the initial 7r-pulse is broadened in fre- 
quency, and therefore the number of excited ions gener- 
ating amplified spontaneous emission in the original fre- 
quency window is reduced. Secondly, there are ground- 
state ions Stark-shifted into the detection frequency win- 
dow which increases the absorption. Assuming a perfect 
TT-pulse with a 1 MHz width, and an optical depth of 1, 
the expected reduction of noise due to these two effects 
is only ~5. However, in practice the 7r-pulse is far from 
perfect and most of the ions in the ensemble see a pulse 
area greater or less than tt (due largely to the Gaussian 
spatial profile of the pulse intensity) . This and structure 
on the inhomogeneous line results in a large free induc- 
tion decay (FID). The third and most significant effect 
of the applied electric field is the reduction of this FID. 



6 




-10 10 20 30 40 50 60 70 



time (us) 

FIG. 7: (Color online) Noise level at the heterodyne beat 
frequency between applied pulses and the local oscillator 
(6 MHz) without any electric field applied (solid green), and 
with the field applied (thin, dashed black). Also shown is the 
shot noise level (dot-dashed red), and the electrode voltage 
profile (thick, dashed blue - with the scale on the right). The 
TT-pulses occur at and 42.7 /is, and the noise level in the 
period between the two pulses is seen to be dramatically re- 
duced when the electric field is applied. If the vr-pulses were 
perfect, the noise following the latter pulse would also be close 
to the shot noise level. For this data series, aL = 0.15. 
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FIG. 8: (Color online) The fraction of atoms in the excited 
state shortly after two 7r-pulses. The medium is initially in 
the ground state then two Gaussian shaped 7r-pulses are se- 
quentially applied to the medium. The pulses were separated 
by 20 times the FWHM of their amplitudes. The plot is gen- 
erated from the solution of the one dimensional semi-classical 
Maxwell-Bloch equations. 



VI. DISCUSSION AND CLASSICAL 
SIMULATION 

There are two areas where these initial experiments 
need to be improved before the proposed memory can 
fulfil its theoretical promises. The demonstrated memory 
is of low efficiency and noisy. The noise measurements, 
Fig. [71 show that the broadening field is successful in 
decouphng excited atoms from the output field while they 
are inverted, but there is excess noise in the time region 
where the light would be recalled from the memory. 



The predominant reason for this noise is the imperfect 
TT-pulses resulting in random FIDs due to structure in 
the inhomogeneous line. In systems that have long lived 
spectral holes like Pr^"'":Y2Si05 this structure is hard to 
avoid. A four-level echo protocol provides a way around 
this problem [l^. Experiments involving Tm'^+iYAG, 
which has short lived holes, have achieved 7r-pulses with- 
out measurable random FIDs 23] . 

Even in the absence of structure on the inhomogeneous 
line, this HYPER quantum memory scheme is still reliant 
on doing 7r-pulses well, but less sensitively. In the output 
region the solution for the time domain is an expression 
like Eq. [31 but with the subscript 1 replaced by 3, for 
region 3. Here a^{Q,t) represents the light incident on 
the crystal during the recall of the light (usually the vac- 
uum). The term that includes 2?3(z,t) represents that 
light emitted from the sample due to any excitation in 
the crystal after the second 7r-pulse. In the case of per- 
fect TT-pulses and vacuum input, 'D^[z,t) will represent 
ground state atoms. In the case of imperfect Tr-pulses 
the atoms will have some unwanted excitation resulting 
in noise. The ability to appl y a ccurate Tr-pulses has al- 
ready received careful study [H, . Semi-classical 1-D 
numerical calculations, presented in Fig. [SI show that in 
the case of no transverse variation in the optical inten- 
sity, atomic population very close to the ground can be 
achieved over a significant fraction of the bandwidth of 
the pulses. This in turn means that the amount of noise 
in the output would be reduced to much lower than one 
photon per spatio-temporal mode. The fact that there 
will be spontaneous emission while the atoms are excited 
will be another noise process, but like the effect of im- 
perfect TT-pulses this noise would not be phase matched 
with the echo. As long as the population is predomi- 
nantly in the ground state when the echo is retrieved, 
the amount of noise will be much less than one photon 
per spatio-temporal mode. 



VII. CONCLUSION 

We propose and present initial experiments for a new 
quantum memory protocol. The proposal uses strong 
rephasing pulses rather than structure in the inhomoge- 
neous line made with spectral holeburning. Noise associ- 
ated with echoes forming in the excited state is avoided 
by using electric field gradients to generate artificial inho- 
mogeneous broadening while the media is in the excited 
state. We present a demonstration of these echoes as 
well as initial noise measurements. In these experiments 
the echoes were noisy due to random FIDs arising from 
structure on the inhomogeneous line rather than ampli- 
fied spontaneous emission from the excited state atoms. 
In order to make the memory quiet these random FIDs 
must be avoided. Furthermore to ensure the atoms form 
the echo while close to the ground state accurate TT-pulses 
will be required. Both of these steps should be techni- 
cally feasible. 
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Note added: Recently, a new proposal has been de- 
veloped that uses some of these ideas [1^. This new 
proposal is attractive in that it removes the need for aux- 
iliary broadening fields. 
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Appendix: Backward Retrieval 

Here we show that similar to other proposals [l^, EB, 
[26j |. in the limit of large optical depth the memory has 
perfect efficiency when the echo is retrieved from the sam- 
ple in the backward direction. 

Here we consider a sample that lies from z = to 
z = L. Equations Hal and llbl are equations of motion for 
the atomic (-Di) and optical (ai) modes propagating in 
the forward direction, toward larger values of z. For clar- 
ity the labels Df and a/ will be used for these quantities. 
The equations of motion for the counter-propagating, 
'backward' fields are the same as the forward case ex- 
cept for a minus sign in Eq. llbl 

^ab{z,t) = -— dADt (A.l) 

The treatment in this appendix is semi-classical. This 
simplifies the expressions greatly because the amplitudes 
of modes with no excitation can be set to zero, and then 
ignored. A quantum treatment will yield the same re- 
sults, but with the appropriate addition of vacuum modes 
to preserve the commutation relations in the presence of 
loss. 

We consider an input field a/(0, t) — ain(i) that is zero 
for t > 0, and atomic modes that are initially in the 



ground state. In this case, we can see from Eqs. [l]and[3] 
that the forward propagating atomic modes at t — are 
given by 



Df{z,t^O,A) = i / af{0,t')e-'^e-'^'' dt' 

-/ — OO 

= ie"T^ain(w = A). (A. 2) 

Here ain(w) is the Fourier transform of ain(t). The upper 
limit on the integral can be changed to infinity because 
the input field is zero for t > 0. 

We reduce the dynamics during the rephasing period, 
the time period between when the broadening field is 
first turned on and the second rephasing pulse, to an 
instantaneous operation at t = 0. For the case of counter- 
propagating TT-pulses this has the effect 

Db{z, t = 0, A) = Df{z, t = 0, A)e-^^*'' (A.3) 

The time delay, td, is the difference between how long 
the atoms are in the excited state and the ground state, 
during the rephasing period. Analogous to Eq.[2]we have 
for i > 



ab{z,uj) = ai,{L,uj)e 2 —ia / dz'e "^^^ Db{z' ,0,uj) 

J z 

(A.4) 

The first term on the right-hand side can be ignored 
because no light is incident on the z = L end of the 
crystal during the retrieval process. The integral can 
easily be computed, resulting in a value for aout(i^) = 
af,(0,w) of 

-ai„(L.)(l-e-"^)e-^"*'^ (A.5) 

or alternatively 

aoutit) = -ainit - td) (1 - e-"^) , (A.6) 

leading to an efficiency that varies with optical depth 
as (1 — e~"^)^. This is equal to 1 for large aL. 
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